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Abstract 

We provide a method to study hadronic matter at finite density in the context of 
the Sakai-Sugimoto model. We introduce the baryon chemical potential through the 
external U(l) v in the induced (DBI plus CS) action on the D8-probe-brane, where 
baryons are skyrmions. Vector dominance is manifest at finite density. We derive the 
baryon density effect on the energy density, the dispersion relations of pion and vector 
mesons at large N c . The energy density asymptotes to a constant at large density 
suggesting that dense matter at large N c freezes, with the pion velocity dropping to 
zero. Holographic dense matter enforces exactly the tenets of vector dominance, and 
screens efficiently vector mesons. At the freezing point the p — tttt coupling vanishes 
with a finite rho mass of about 20% its vacuum value. 



1 Introduction 



Recently there has been much interest [HE] in the AdS/CFT approach [3] to study the non- 
perturbative aspects of gauge theories at large N c and strong coupling with even applications 
to heavy ion collision experiments such as jet quenching [3]. On the other hand, dense QCD 
is of much relevance to compact stars and current heavy ion collisions. This subject has been 
intensely studied in the past decade following the observation that at asymptotic densities 
QCD matter may turn to a color superconductor due to asymptotic freedom [5j. At very 
large N c there is evidence that the Overhauser effect takes over with the formation of density 
waves [Bj. Indeed, at weak coupling the high degeneracy of the Fermi surface causes quark- 
quark (antiquark-antiquark) pairing of the BCS kind with different color-flavor arrangements, 
while at large number of colors quark-antiquark pairing of the Overhauser kind is favored. 
In strong coupling, first principle calculations are elusive and the two pairings may compete. 

To address the issue of strongly coupled QCD at large N c in dense matter we explore in 
this paper the holographic principle. In recent years, the gauge/string duality has provided 
a framework for addressing a number of problems in strong coupling where little is known 
from first principles in the continuum. To do so, we consider a model for holographic QCD 
recently discussed by Sakai and Sugimoto. In the limit Nf <C N c , chiral symmetry in 
QCD is generated by immersing Nf D8-D8 into a D4 background in 5 dimensions where 
supersymmetry is broken by compactification (Kaluza-Klein mechanism with M^k scale). 
The induced DBI action on D8 yields a 4 dimensional effective theory of pions and infinitly 
many vector mesons where the M^k scale plays the role of an upper cutoff (the analogue of 
the chiral scale Att f n ). At large N c baryons are Skyrmions. 

In this paper, we extend the analysis by Sakai and Sugimoto [2] to finite baryon density. 
Although there have been many papers [7] on finite chemical potential in the holographic 
approach, all involve the R-charge or isospin instead of fermion number. The difficulty for 
the fermion number is that the U(l) charge for R-symmetry is already dual to the U(l) 
charge of AdS Reisner-Nordstrom Black hole charge. There is no additional U(l) in bulk to 
modify the geometry of the AdS black hole. In this paper we follow the more traditional 
route of introducing the chemical potential via the nonperturbative induced DBI plus CS 
action through the external U(l)y source in Sakai and Sugimoto chiral model. Specifically, 
we introduce the baryon chemical potential as V = —i^i^l^c where is the external vector 
field in the induced DBI plus CS action. While Vo drops in the DBI part it contributes to 
the CS part through the WZW as it should. A redefinition of the vector fields in the Sakai- 
Sugimoto model necessary to enforce vector dominance, causes a reshuffling of V from the 
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CS to the DBI action as we detail below, with important phenomenological consequences. 

In section 2 we review the emergence of the induced effective action on the probe D8 
brane from the 5 dimensional D4-brane background, and we detail the contributions of the 
DBI action to fourth order in the vector fields, as well as the Chern-Simons part. The 
expanded effective action enforces total vector dominance. In section 3 we introduce the 
baryon chemical potential and details the nature of the effective expansion at finite density. 
In section 4, we extract the baryon number density, the pressure and the energy density 
of holographic dense matter. In section 5, we discuss the pion and leading vector meson 
parameters and interactions in holographic dense matter. Our conclusions and suggestions 
are summarized in section 6. 

2 Duality and Branes 

The generalized holographic principle states that at large N c maximally supersymmetric and 
conformal gauge theories in flat D-dimensions are dual to Superstring/Mtheory on pertinent 
AdS spaces, i.e. for D = 4 it is AdS 5 x S 5 and for D = 6 it is AdS 7 x S 4 . 

To model QCD in D = 4 dimensions, we start from D = 6 superconformal theory, and 
compactify it twice. First through S* 1 with radius Ri leading to D = 5 SUSY theory. Second 
through another S* 1 of radius R 2 ^> Ri with fermion antiperiodic boundary conditions to 
break SUSY. Below the cutoff scale M KK = I/R2, the D=5 broken SUSY theory behaves as 
a D=4 gauge theory. 

2.1 D4/D8 Branes 

Doubly compactified D=6 is the boundary of AdSj in which the M-theory reduces to type IIA 
string theory. As a result M5 branes wrapping around the first S 1 with radius R 1 transmute 
to D4-branes. iV c copies of these branes yield near extremal black-hole background. 

The metric, dilaton and the RR three-form field C 3 in a D4-brane background are given 

by 

ds 2 = g) 3/2 (rj^dx" + f(U)dr>) + (^) 3/2 (j^L + , 

^ = 9s{^)'\ F 4 = dC 3 = 2 -^e 4 , f(U) = l-^. (1) 
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Here /i — 0, 1,2,3 and r is the compact variable on S 1 . U > Ukk is the radial coordinate 
along 56789, e 4 the 4-form and V 4 = 8tt 2 /3 the volume of a unit S 4 surrounding the D-4 
brane. R 3 = -ng s N c l 3 , where g s and l s are the string coupling and length respectively. Again, 
this background represents N c D-4 branes wrapped on S 1 . 

Let us consider Nf Z)8-probe-branes in this background, which may be described by 
t(U). We choose a specific configuration such that r is constant(r = 4 m )■ ^ corresponds 
to the maximal asymptotic separation between D8 and D8. with a new variable z, instead 
of U, defined by the relation, 

U={U 3 KK + U KK z 2 )^=U z (2) 

the induced metric on D8 is 

ds 2 m = gMN<ix M dx N 

= (l) 3/2 '^^(6 3/2 ^ 2+ (l) 3/2f/M - (3) 



where 



M{N) = {/x(z/)(0,l,2,3),z(4),a(5,6,7,8)} (4) 



2.2 DBI and CS Action on D8 

Consider the U{Nf) gauge field A M on the probe D8-brane configuration. The effective 
action is 9-dimensional and reads 



Sds — Sdbi + Scs 



Sdbi = -T J d 9 xtr (e~V- det(g MN + W^kj) (5) 
Scs = C 3 trF 3 = -L-[ F 4 u 5 (A), (6) 



a 

where = g s (^r) 4 , F 4 = <iC 3 is the RR 4-form field strength and lj 5 (A) is the Chern- 
Simons 5-form, 



u 5 (A)=tr[AF 2 - 1 -A 3 F+±A 5 ) , (7) 



which satisfies duj^ = trF 3 . 
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Assuming that A a = and and A z are independent of the coordinates on the S A , the 
induced DBI action becomes 5-dimensional 



Sdbi = -T J d 9 xtr (e^y/- det(g MN + 2na'F MN ] 



T Id xdz U, 



m 9 77 

tr J 1 + (27ra') 2 -—r ] ^r ] P-F^F ua + ^a') 2 --^^ F^ Z F VZ + [F 4 ] + [F^ 
lUX 4 (Jkk 



where T = . [F 4 ] and [F 5 ] are short for structures of the type 



(9) 



[F 4 ]-F^u + ^ u (10) 



The Chern-Simons term is 



S cl = ^ f , (12) 



with the normalization ^ / 54 F4 = iV c , and M 4 x R is the five- dimensional plane parameter- 
ized by x^ and 0. 

2.3 Effective Action in 4-dimension 

In this part we review how the 5-dimendional DBI action yields the 4-dimensional effective 
action. Essentially, the 5-dimensional induced action with compact S 1 is 4-dimensional for all 
excitations with wavelengths larger than the compactification radius of the order of 1/Mkk- 
Throughout and for zero baryon density, our discussion parallels the original discussion by 
Sakai and Sugimoto [2]. We quote it for notation and completeness, and refer to their work 
for further details. 

The leading terms in the 1/ A ~ 1/ Mkk expansion of the DBI action is 



S™ = K I d A xdZ tr 



-K~ 1/3 F 2 +KF 2 

2 jW 1 fjLZ 



(13) 
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where 

« = f (2W) 2 i? 3 = , Z = i^l + Z 2 . (14) 

1U07T° UkK 

In order to extract four-dimensional meson fields out of the five dimensional gauge field, 
we expand the gauge field as 



A,(x^z) = , (15) 

71=1 

OO 

A z (x»,z) = ^ (O H^o(z)+Xy n V)0n(z) , (16) 

n=l 

We choose the functions ip n {z) to be the eigenfunctions satisfying the self-adjoint differential 
equation 



-K l ^d z {Kd z ^ n ) = \ n ^ n (17) 

where A n is the eigenvalue and {ip n } is a complete set. the normalization is fixed by the 
Laplacian in (flj 



k J dZ K 1/3 ip n ip m = 5 nm . 

This implies 



k J dZ Kdz^ndz4> m = A n 5 nm (19) 
The 4> n (Z) are chosen such that 

4> n {Z) = (m n U KK )- l d z ij n (Z) (n > 1) 

MZ) = r-Jrj K (20) 



satisfying the normalization condition: 

9 - T{2ira') 2 U* K J dZK<f> m <j) n = 5 mn , (21) 



■'mi fn) 



which is compatible with (fT7|) and (TT8j) . 



5 



Inserting (fl~5]) and (|T6|) into (|T3|) and using the orthonomality of ^ n and n , yield 



S 



DBI 



d x tr 



(0)\2 



oo , 
n=l \ 



(n) 
A* 



+ (interaction terms) . 



(22) 



In the expansion ( [151) and (TIBl) . we have implicitly assumed that the gauge fields are zero 
asymptotically, i.e. Am{x^,z) -> as z -> ±oo. The residual gauge transformation that 
does not break this condition is obtained by a gauge function g(x M , z) that asymptotes a 
constant g(x M , z) — > <7± at z ± oo. We interpret (g + , gJ) as elements of the chiral symmetry 
group U(Nf)i x U(Nf)n in QCD with Nf massless flavors. 

2.3.1 External photons 

By weakly gauging the U(Nf)i x U(Nf)n chiral symmetry, we may introduce the external 
gauge fields (Al m , Ar m ). For the interaction between mesons and photon A° m , we may choose 



A-Lfj, = A RfM = eQA 



em 



(23) 



where e is the electromagnetic coupling constant and Q is the electric matrix-valued charge 



Q 



( 



\ 



-l 



/ 



(24) 



for Nf = 3. To insert the external-source gauge fields, we impose the asymptotic values of 
the gauge field on the D8-probe-brane as 



lim A^af, z) = A Ltl (x") , lim A^x^z) = A Rfl (x>") 

z — ►+oo 



z— > — oo 



(25) 



This is implemented by modifying the mode expansion (II 5p as 



A^(x^, z) 



A Ltl (x^+(z) +^Si B) MVv. 



n=l 



where ip±{z) are defined as 

1 2 

V'iW = oi 1 ± Vw) , ^o(«) = - arctan z , 

Z 7T 

which are the non-normalizable zero modes of (flTI) satisfying d z ip±(z) oc 0o(^)- 



(26) 



(27) 



2.3.2 A z = gauge 

The A z = gauge can be achieved by applying the gauge transformation Am — > gAug^ 1 + 
gdhig' 1 with the gauge function 

^(xV) = Pexp j-jf dz' A z (x»,z')^ . (28) 



The asymptotic values of (1251) change to 

lim A,(x», z) = A%(x") , hm A,(x», z) = A%{x») , (29) 

2— >+OC r Z^ — OO ^ 

and 

= e+M^MtfM + ^+M^tfM. (30) 

4-(x^) = e-M^wc'w+e-W^- 1 ^)- (31) 

where £±(x M ) = lim^-too g(a; M , 2). The gauge field in the A z = gauge can be expanded as 

00 

A^x^z) = A%(x^M^)+A%(x^(z)+J2B { ;\x^ n ^) ■ (32) 

n=l 

The residual gauge symmetry in the v4 z = gauge is given by the z-independent gauge 
transformation. The residual gauge symmetry h(x^) G U(Nf) and the weakly gauged chiral 
symmetry (g + (x M ), g_(x M )) G U(Nf) L x U(Nf) R act on these fields as 

-4l m -> 9+Al/a9+^ + 9+d»g+ l , (33) 
-> 3-Zl^Z 1 + tf-^Z 1 , (34) 

e± -> ^Ci^i 1 , (35) 

5( n) -> hB^h- 1 . (36) 

Here £±(3^) are interpreted as the U(Nf) valued fields £,l,r{x^) which carry the pion degrees 
of freedom in the hidden local symmetry approach. Actually the transformation property 
(j35p is the same as that for ^l,r{^) if we interpret h{x^) G U(Nf) as the hidden local 
symmetry. They are related to the U(Nf) valued pion field U(x^) in the chiral Lagrangian 
by 

ezV)^(^) = U{s?) = e 2m ^ )/f * . (37) 

The pion field IT(x M ) is identical to (p^\x^) in ( TlBl) in leading order. 

1 Here the pion field 11(2;^) is a Hermitian matrix, while (p^{x^) and the vector meson fields BjrHx^) 
are anti-Hcrmitian. 
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If we choose the gauge 



= 1, e; 1 ^) = U(x») = e 2iU ^ )/f * (38) 

then the gauge field is 

A^x^z) = U-\x»)A Lll {x»)U{x^ + {z)+A R ^)ij„{z) 

+ U-\x»)d,U(x»)i, + (z) + Y,B { ;\x^ n (z) (39) 

n>l 

An alternative gauge is also 

^(x") = = e m{x ^ /f * . (40) 

In this gauge, the gauge potential in (1321) can be expanded up to quadratic order in fields as 

K = (^ + ^[n,9 M n]-^[n,^]) + ^ + ^ M n-^-[n,v M ]^o + 

oo oo 

+ ^<^ 2n -i + ^a^ 2n + --- , (41) 



n=l n=l 



with 

V,= ^(A LlM + A R ,) , A ^ l -(A L ,-A R ,) , = St*- 1 ) , < = (42) 
2.3.3 Skyrme action and WZW term 

It is interesting to note with Sakai and Sugimoto [2] that the induced effective action ex- 
panded to fourth order in the gradients yields exactly a Skyrme-like action 



where the pion decay constant /„. and the dimensionless parameter es are fixed 



, (43) 



- ^=ih XN - (44) 

ef = k fdzK- 1/3 (l-^) 2 . (45) 
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Also in the A z = gauge, the Chern-Simons term is generic 



at r 

^ 47r J M 4 



+ - > 9d9 ~ I)5 1 ■ (46) 



where g is the gauge function given in ([28"]) and the 4-form «4 is 

a 4 (^,A) = -^tr^(ArfA + rfAA + A 3 )- VaVA-^ 3 ^ . (47) 

As expected ( l4T)j) yields the WZW term in QCD with gauged pion fields. Indeed, if we omit 
the vector meson fields Bj?\ we have 



S, 



cs 



X " ' Z ~^tJ tr( ^" )5 ' (48) 
M 4 2A0tt z J M 4 xR 



„«=a"=0 487T 2 J M 4 240tT 2 



where 



Z = tv[(A R dA R + dA R A R + A^U^AlU + U- l dU) -p.c] + 



+ tr[dA R dU- 1 A L U - p.c] + trfAj^dtr 1 ?/) 3 - p.c] + 

+^ trader 1 */) 2 - p.c] + trfC/A^tr^dC/dC/" 1 - p.c] - 

- trlARdU-'UARU-'ALU - p.c] + i tr^C/" 1 ^) 2 ] . (49) 
Here "p.c." represents the terms obtained by exchanging Al Ar, U <->• /7 -1 . 

2.3.4 Vector meson dominance 

The induced DBI action f fT3|) carries exact vector dominance. Indeed, by diagonalizing the 
kinetic terms of the vector meson fields, vector meson dominance emerges from the underlying 
DBI action naturally. To see that, we recall that the gauge fields are 



A„ = V^ + ^o + X^^-i + X^^, (50) 

71=1 71=1 

A z = IT 0o, (51) 
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and involve both the external vector sources V, A and the dynamical vector fields V, a. To 
diagonalize the kinetic terms of the dynamical vector fields require the introduction of the 
physical tilde vector fields v, a which are 



in terms of which the gauge fields now read 



n=l 



n=l 



with 



J 2n ■ 



n=l 

J d,Z K~ 1/3 tp2n-l , dAa™ = K j dz K^^Zn^O 



n=l 



The DBI action (1131) in terms of the physical vector fields to fourth order is 



(52) 
(53) 



(54) 



(55) 



K J dz tr 
1 



2 ^ 



tr 



l ^((^) 2 + (^) 2 ) + 

+ {d^a* - cU£)([^, a 9 1 - a?»]) g m + 

+ \ K, <} ^1 9v^ vPvq + l - [a™, al) [a***, a qv ] + 

+ ([<,^][a PM ,a^] + [t^a£][^a<H - p?,fl?][5 m ' l a*"]) g mPaq 



(56) 



Note that all the couplings between the external gauge fields (A L , A R ) and the vector meson 
fields (v n ,a n ) vanish in the first term of the effective action ([13]) . The second term in the 
effective action (TT31) 
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njdztx [KFl] 



tr 



n 2 vn (v; - a Vv nV^ 2 + m 2 an (a; - aj^Atf + (id„U + UA p ) 2 + 



z a n am [n, a;} [n, a"*] - c„ [n, ^] [n, 



(57) 



The mesons couple to the external gauge fields only through v n — > V and a n — > A transitions 
in f7]). 

Finally, inserting (1501) and (1511) into the Chern-Simons action gives 



S 



cs 

D8 



N c % 



tr 



UdB n dB m c nm + 



4vr 2 U j M 4 

+ II (dB m B n B p + B m B n dB p ) c mnp + UB m B n B p B q c mnpq 



+ 



+ 



24tt 2 



tr 



A/ 4 



B m B n dB p d mn \ p - -B m B n B p B q d 



mn.p| q 



(5? 



where i? 



2n-l 



f n , = a n It shows complete vector meson dominace in the WZW term. 



The terms with two or more pion fields are absent. 



3 Induced action at finite density 

For simplicity we work in Euclidean space, in which case the DBI action reads 



S 



E 

DBI 



T I d 4 x E dzU 2 



tr W 1 + (W)«|^F*F* + (2W)^^F£F£ + [(F*)*] + [(F*)«] 



; KK 



(59) 



following prescription. 



o 



X 

S E 



-IT , 

-iS . 



(60) 
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From now on we omit the superscript E. From the arguments to follow, we will see that due 
to vector dominance no matter term is generated through the Chern-Simons part. In a way 
this is expected; the anomaly in the Chern-Simons form is an ultraviolet effect that after 
regularization shows up in the infrared. Since matter dwells mostly in the infrared, the scale 
decoupling insures the insensitivity of the Chern-Simons terms therefore of the anomaly. 
This result is usually referred to as the non-renormalization theorem. 

The density will be introduced through the boundary and via the external vector field V 
through 

V M = -i/^Sfjo l Nf xN f , A^ = 0. (61) 

The master gauge field now reads from (|50|) 



n=l n=l n=l 

For calculational convenience below, we split Vq in two parts: a constant U(l) part (— Wq C ) 
and the rest (Vq G ) , 



= -Wq - ia n fi = -Wq C + v% c - ia n \i 



where —Wq C and Vq C are the constant U(l) part and the rest part of respectively. Thus 

oo oo oo 

n=l n=l n=l 

oo oo 

A = 5>r^2n-l + X>"^2n > ( 63 ) 
n=l n=l 

Defining, 



j^2n-l 
X 2n 



B 



2n-l 




B 



2n 







v, 



oc 



B 



2n-l 



(64) 
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leads 



71=1 

A z (x»,z) = -iU(x)(j>o(z) . (65) 
The corresponding field strengths are 



F, v = J2( d ^-duB;)^ n + [B^K]^ m , (66) 

n=l n,m=l 

oo oo 

f z , = id^n 0o + + b;(x») + ^[n, b;]^ . (67) 

n=l n=l 

Notice that iX n 5 fJi0 does not contribute to F^ u . The additional terms due to iX n S IJi0 come 
from F^ z . They contribute to the induced action through 



F^F^ = {d^-d v B^){d,B™-d v B™)^ m 

+ {[5;, B™\ , (d,Bl - d u Bl) } M m ik , (68) 

f^ z = -(d^n)<pl + B;B^ n1 p m -x n x m ^ m 

+ (i{d,n , b;} - 2(d u)x n ) o ^ n + 2ix n B™i>J m 

+ (-2X n [n , b™} + i{b; , [n , £>]}) ^ m <p + {^n , [b; , n] } <\>^ n , (69) 

[F 4 ] = (2na') 4 ljj£^ {iWWo + (iX n + B^} l -F t] F t] {*(d o n)0 o + {iX m + F m )^ m } 

[F 5 ] ~ (2W) 5 X"X m ^ m i(F 12 F 2 3F 31 -F 31 F 12 F 2 3) (70) 

For notational simplicity we have omitted X^Lr 

The induced DBI action at finite density can be separated into field independent and 
dependent parts Pq and P\ respectively 



Po = i-( 2TO ') 2 ^rrw m 

4 L/rk 

Px = (2W) 2 -^F^ + (2W) 2 ^(F M2 F )1Z + X"X m ^ m ) + [(F) 4 ] + [(F) 5 ](71) 
so that 
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with 



S 



DBI 



T / <Fxdz Ut tr J P + P 1 



T / d 4 xdz Ut tr 



Po + r 



S 1 (X n )+S 2 [U,v, ] X n ] + 



1 P l 1 Pp 



2v^ 8^ 



+ 



Si(X n ) 



T / d^xdzU; tr W 1 



T / cfxdz P 2 tr 



1 A 
2 



4C/ KK 

1 Pf 



+ 



(72) 
(73) 



We note that the dependence on the chemical potential /i in Si is of the the type 



\J\ — j^p?. There is no term linear in //(except through the Chern-Simons term which 
is odd under t — > —t). The reason is that the QCD partition function is even under fi — > —\i 
since the matter spectrum is symmetric around zero quark virtuality. We believe that this 
behavior maybe derived geometrically from a change in the underlying metric, but we are 
unable to show it. 

The chemical potential /x as defined in (j50|) with ([ETC]) (i.e. V M = —i ljVfxJV,, -A^ = 0) 
yields P^ = 0. Thus there is no contribution to the DBI action. The only contribution 
stems from the Chern-Simons action as can be checked explicitly. However the definition 
(I50|) is at odd with VMD as explained in detail in [2J. In other words the physical vector 
field are defined in (I52|) and not t>™ in (150]) . So (150]) has to be substituted by (R?2"|) . In this 
case Pq 2 7^ and \x contributes to the DBI action but not to the Chern-Simons action. 



4 Free energy of dense matter 

On the boundary the induced DBI action describes a dense hadronic system at finite baryon 
density whereby the baryons are Skyrmions in large N c . The free energy Q is given as 



Q = -TlnZ 



G 



— Tin \N' J D[^]D[n]e- s ^ 
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-kTln U'j 2)[^]D[n]e-^( xl ) +52 [ n ^ ;X "]) 
TSi - Tin (j 2)[^]D[n]e- 52[n '^ ;X " ]) 



(74) 

where is the grand partition function and N' is an irrelevant matter independent constant. 
The Gibbs relation 

Vl = U -TS - fiN = -PV (75) 

implies that Q(T = 0) = U—fiN is the free energy at zero temperature and finite baryon den- 
sity. Ignoring the meson zero-point contributions (classical limit), the leading contribution 
to the free energy per unit volume (minus the pressure) is 

- Q ( T = °) gi(* w )-Si(0) 

e{x ] = ~ — m — 



T I dzU 2 z ti (2n a ')^^-X^ijJ m - 1 



rn 



W .-S^-.l, (76) 



27 3 tt 5 ./ \ V N C M 2 KK ^ 



where AT = 1 + Z\ X = g 2 YM N c , and ^ n = ^ n {Z) = ^ T(2TiaYR^ n (Z) . 
In terms of the meson physical parameters fixed at zero density^ 

2 A r f 2 277T 4 



M 2 KK = ^ k = 1,2,3, 
the free energy density becomes 



with 



iV c m 10 J \ V U lJirO 



iV c ' 

Cj = A l7 r 7 , Dx = A??r 7 . (77) 



2 From this section the subscript in the physical quantity, like f^o, means the value at zero density. 
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The sum over the vector meson-species m, n is subsumed. ( 1771) receives corrections from 
the quantum loops with mesonic insertions. In a way (!77|) captures the deformation of the 
underlying D8-probe-brane in the presence of a fixed External Vo to leading order. The 
square root shows that the response is not analytic, with a subtle branch point singularity 
to be physically interpreted below. 

In the limit N c ^> 1 but fixed baryon density [is = ^N c (low density limit), 



im 1Q N c ( n fiB ~ n \ / 'm^B 



Ks w a c a n — - v" r ) \ a m — - v™ ^ 2n _i# 2m -i 
-L(^/N~ c l- _ 5» )2 _ J_ ( 78 ) 



The free energy density can be Tayor expanded as 



N f M KK ( n n^B ~n \ ( n m^B ~ m I / , y , - ,t. ,t, 

a — V 0C ) \ a Tt V 0C / d ^ K ^2n-1^2m-l 



Nf A o , /T7- 1 ™ ,o -Vr 



/X) m (2n-l)o(V^" 



"'(*.-l)0VV-c^-«0OJ 
n=l 

which is quadratic in /i# and of order Nf/N c which is small. 

In the limit N c ^> 1 but fixed quark density \i = N c / fig (high density limit), 



.1 mj N, 



[a fx 



J oc) 



*2n-l* 



2m- 1 



(79) 



so that 



/V 00 

~ -^E^-Do^) 2 -^/^ ( Xn<<1 ) (80) 

71=1 

which is of order NfN c . 

In Fig. [T] we show the free energy versus X 1 /m p o for N c = 3, Nf = 2 and m p o = 776 MeV, 
f n0 = 93 MeV. Since 

2 

777. 

fc0 ^A fe = 0.67, 2.88, 6.6, 11.8, ■■• (81) 
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The baryon density n B follows from the free energy through 



n B 



n b i on 



N,Ml K ( dZK - ff| 

J . A _ Tfh<o N h 



f dZ , c ^ ~ JV° fixed a B 

° ^ (82) 

fdZ ~N C fixed /i 

where again only the n — m — 1 contribution was retained. This is a transcendental equation 
that reexpresses the shifted chemical potential X in terms of the physical baryon density 
n B - For large N c and fixed fi B it shows that n B ~ N® which we interpret as the low density 
phase. For large N c and fixed \x = fi B /N c it shows that n B ~ N c which we interpret as 
the high density phase. In Fig. [2] we show the numerical solution to ([8*2]) for n B /no versus 
X/rripo, where tiq = 0.17fm 3 is nuclear matter density. Asymptotically large densities are 
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attained for a critical X/m p o m 0.05 \fTJ\ f^ /N c m^ . For N c = 3 and Nf = 2 we have 
X/m p Q ~ 0.111 and n B /n ~ 43. 

The existence of a limiting chemical potential reflects on the fact that the underlying 
D8-probe-brane becomes unstable if the external V becomes very large, as the square root 
in the ground state energy developes an imaginary part in. 




0.12 



Figure 2: ^ vs — 

° no rripo 



The pressure is P = — e which is 



In Fig. [3] we show how the pressure changes with baryonic density ub/uq for N c = 3, 
Nf = 2, m p0 = 776 MeV and f w0 = 93 MeV. The pressure is dominated by repulsive two- 
body Skyrmion-Skyrmion like interactions at low density (the term of order n 2 in (153]) ). At 
high density attractive three-body Skyrmion-Skyrmion-Skyrmion like interactions appear and 
cause the pressure to saturate (the term of order n 3 in (1531 ). The two-body interaction are 
100 times stronger throughout the density range explored. The specific fit is 
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P = (1.17 10 7 ) n 2 - (2.14 10 5 ) n 3 



(83) 



If we recall that in matter the kinetic energy K = (dP/dn) /2, then the slope in Fig. [3] is a 
measure of the kinetic energy. The holographic principle suggests that hadronic matter at 
large N c freezes at % ~ 30 no- 




Figure 3: Pressure vs — 



The energy density is defined as 



' v 3 v 3 

= e + fi B n B = e + —^{vl c + X x )n B 
a 

= e + ^MnB + ^fe (84) 

and is sensitive to the gauge choice of the dynamical U(l) v vector meson (omega). In the 
presence of U{l)v mesons, v™, the energy density is dependent on the constant part of the 
U(l)v meson field, i.e. Vq C . This is reminicent of the gauge sensitivity of the EM potential. 
In the gauge where Vq C = 0, the energy density is shown in Fig. HI The repulsive three-body 
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interactions soften the rise in the energy density. Note that the rise becomes almost linear 
after n B /n ~ 30. The slope is sensitive to the choice of the gauge v\ c = 0. Since the 
pressure curve shows total freezing at these large densities, we expect that hadronic matter 
crystallizes. The slope is expected to be negative and a measure of Madelung constant for 
the specific crystal symmetry. 




Figure 4: Internal energy density vs — 



5 Mesons in Holographic Dense Matter 

The DBI effective action in matter fixes completly the meson dispersion laws and interactions. 
Indeed, in matter S^ff^X™] is 



S 2 [U,v^X n ] = / d 4 xdZtr 



QTU KK U 2 Z P 1 - - A Q*TU KK U 2 z Pl + 



15) 



where 



1 



1 - AKr^j 
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Only the term n = m = 1 in the summation over the vector meson species was retained (See 
(15Tj) ). In this section we treat as ant i- Hermit iaon matrices with generators t a normalized 
as tr(t a t b ) = |<5 a fc. The effective action is rotated back to Minkowski space. Throughout, 
only terms up to the third order in the fields are retained , i.e.O((IT, t^) 4 ). The leading terms 
to this order in [F 5 ] vanishes because of the cyclicity of the trace. 
The expansion can be simplified by using the notations 



v := v. 



i 



X :=X l 



*7) 



so that 



S 2 [U,v fJr ;ii] 



d x 



a£ 2 tr (<9 n<9 n) - a£ 2 tr (^n^II) 
1 



- a v 2 tr(d Vi - diV ) + -a v2 tr(<9;^ - djVi) 

- m v tr v + m v tr v { 

- a% 3 tr ([v , Vi](d Vi - <9^ )) + tr ([vi, Vj)(diVj - djV { 

- a T vlP tr (d Tl[v , II]) + a s vlP tr ($1%, II] ) 

- a l v-> tr (vo(doVi - <9^ ) 2 ) + a? „ 2 tr (v (diVj - djVi) 2 ^ 

- a-v iP (d Ild Ilv - dottdilLvi) - a VQV 2V (v 2 - vf) 

+ dZKQfavA tr ^(^n)^ + /i^(<9 n) + ^[n] 3 



The last two contributions in (1551) vanish because i^, Q, ^i and 0o are even functions and 
is an odd function. The new two terms related to vq are in the 6th and 7th lines of 
( I88p . All the parameters appearing in (1881) are tabulated below for non-zero density. Their 
respective values at zero density are recorded on the right-most column for comparison. We 
note that Q, A ~ iV c for fixed fi, so that all coefficients are of order N® for fixed fi at large 

Nr.. 
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Coefficients 


nonzero X 


x = o(g = i, a = o) 


T 


i / dzx- 1 Q (l + Ag 2 A'H 2 ) 
i f dZK~ x Q 


i 
i 


T 

a 2 


JdZK-^lQ 
fdZK-^lQ- 1 


i 
i 


2 T 
2^ 


/ c/Z A^ 2 g (i + AQ 2 A3 ^ 2 j 
^fi [dZK&Q 


m 2 




J^fl fdZK-^lQ 
^ [dZK-^tQ- 1 


Zpt ■ 0.45 
^? ■ 0.45 


T 

<n 2 

%n 2 


^= / dZ K~ l ^\Q (l + Ag 2 A3*f j 


^=•0.63 
-t= ■ 63 


T 

a 2 
a s 2 

d V0V 2 


#/^g 3 A-§vi> 2 
|/cizg 3 Al^ 









5.1 Meson Velocities 



For N c = 3, m vo = 776MeV, = 93MeV, we show the pion in Fig. [5] and vector meson 
in Fig. [6] renormalization constants, both for the time and space-components. The resulting 
velocity of the pion and the rho meson in holographic dense matter is shown in Fig. [7J The 
pion velocity 

"'=/i^l (89) 

with fn ,T , the spatial (S) and temporal (T) pion decay constants, approaches zero when 
hadronic matter freezes. The vector meson velocity 

v * = \ji (90) 

is about half. 



5.2 Vector Screening Masses 

The pion is massless throughout as the current quark masses are set to zero in this work. 
The vector meson masses follow by Higgsing, and change in matter both space-like and 
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time-like. Space-like, the vector meson masses are screening masses. The unrenormalized 
masses increase with increasing baryon density as shown in Fig. [HJ The physical screening 
masses (Mr) are obtained by normalizing with the pertinent matter dependent space-like 
wavefunction renormalizations. Both the vector and isovector screening masses increase with 
increasing baryon density as they should. Fig. [9] show the screening masses for longitudinal 
(left) and transverse (right) vectors. They are defined as 



2T / 2 
M 2 = m v / m v0 



T 

a , 2 



ml /m: 



a s 2 

V 



(91) 





Figure 9: Screening masses: Longitudinal (left) and Transverse (right) 
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5.3 Vector Meson Masses 



The time-like vector masses are obtained similarly to the screening masses by instead using 
the time-like wavefunction renormalizations. That is 



Ml = mf /m '° 



(92) 



The dependence of M v on the baryon density is shown in Fig. [TUJ The rho meson mass 
drops by 20% when holographic dense matter freezes. 



[(m^/m^/a^]" 2 











5 10 15 20 25 30 35 40 45 



Figure 10: Time-like Vector mass vs baryon density 



5.4 Vittt and VVV Couplings in Matter 

The Vtttt coupling is modified in matter. The longitudinal and transverse couplings are 
shown in Fig. [11] (left). The longitudinal coupling drops by 40% at the freezing point, while 
the transverse coupling vanishes when dense matter freezes. In Fig. [TT1 (right) VVV coupling 
also decrease in matter. The transverse couplings drop by 90%, while the longitudinal ones 
by 30%. In addition to these standard couplings, new matter dependent couplings emerge. 
Their composition and dependence on the baryon density is shown in Fig. [T2l 

Many of the results presented here in matter bears similarities (and differences) with 
arguments presented at finite temperature using phenomenological models with hidden local 
symmetry [8]. 
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Figure 11: Vittt and VVV couplings. 




Figure 12: New terms 
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6 Conclusions 



We have provided a minimal extension of holographic QCD [2] to dense matter. The induced 
DBI effective action on D8 in the presence of a constant external U(l) v source Vq is the 
effective action we expect from the gravity dual theory from general principles. Indeed, if we 
were to solve for the effect of the boundary Vo on the D4-brane background and the resulting 
DBI action on the D8-probe brane, the answer is the externally gauged DBI induced action 
in the vacuum. The effective action develops an imaginary time at large Vo, a signal that 
the D8-probe-brane craks in the external field. 

Holographic matter describes dense QCD at large N c with baryons as solitons. Its bulk 
pressure asymptotes a constant at large density, signalling total freezing with zero mean 
kinetic energy. Before freezing the matter is dominated by two-body repulsion at low density 
and three-body attraction at intermediate densities. The two-body effects are 100 times 
stronger than the three-body effects. 

In holographic matter the pions stall to almost a stop, while the vector mesons only slow 
down. The vector masses drop by about 20%, while vector screening becomes increasingly 
large. The transverse vector mesons completly decouple from pions at large densities. Many 
of the current results bear similarities with known results from effective models at large N c . 
In a way, they are new as they provide first principle calculations to large N c and strongly 
coupled QCD. 

In the work to follow, we will present results for holographic matter including finite quark 
masses and temperature. 
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